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ABSTRACT 
In the present paper, the mathematical model of single- substrate enzyme kinetics is discussed. We employ the new Homotopy 
perturbation approach to solve the coupled non-linear differential equations containing a non-linear term related to basic Henri- 
Michaelis-Menten rate equation. The analytical expressions for concentration of substrate, enzyme- substrate complex, free enzyme 
and product have been derived for all values of rate constant. A numerical simulation is also reported using Matlab software 
program. Our analytical results are compared with our simulation results. A good agreement is noted between analytical and 
numerical results. 
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1. INTRODUCTION 

The vast majority of chemical transformations inside cells are carried out by proteins called enzymes. Enzymes accelerate the rate of 
chemical reactions (both forward and backward) without being consumed in the process and tend to very selective, with a particular 
enzyme accelerating only a specific reaction. In biochemistry Michaelis-Menten kinetics is one of the simplest and best-known 
models of kinetics. This model contains three rate constant can be derived by temperature jump method or transient state kinetics, 
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but both methods need more complicated techniques and equipment [1]. Wei Chen and Cheng Zhu et al. [2] developed a kinetic 


model for a single-substrate trimolecular enzymatic system. The general kinetics of the trimolecular enzymatic system is more 
complex than the Michaelis-Menten kinetics. The Briggs—Haldane approximation of the irreversible Michaelis-Menten scheme of 
enzyme kinetics is cited in virtually every biochemistry text book and is widely considered the classic example of a quasi-steady-state 
approximation [3]. 

Many modeling approaches rely on ordinary differential equations (ODE) which are based on standard enzyme kinetics. 
Michaelis-Menten and Hill functions are indeed commonly used in dynamical models in systems and synthetic biology because they 
provide the necessary nonlinearity to make the dynamics nontrivial [4]. The most widely used rate expression for single-substrate 
enzyme catalyzed reactions, namely the Michaelis-Menten kinetics is based upon the assumption that enzyme concentration is in 
excess of the substrate in the medium and the rate is mainly limited by the substrate concentration according to saturation kinetics 
[5]. Alberto Maria et al. [6] described that a new asymptotic expansion valid in enzymatic reactions, where the total amount of 
enzyme exceeds greatly the total amount of substrate. In such a case, it is well known that the Michaelis-Menten approximation is 
no longer valid; therefore asymptotic expansion, which improves known results, is a new tool to approximate in a closed form the 
concentrations of the reactants in the presence of an enzyme excess. Time-dependent closed form solutions are derived for the 
three nominal cases of competition: even, slow and fast inhibitors, allowing for the first time the complete characterization of the 
reactions. The time-independent Michaelis-Menten approach is in accurate when a fast inhibitor is present [8-10]. 

Various aspects of the QSSA have been studied from time to time. An early study revealed that the condition is necessary for the 
applicability of the QSSA [11 - 14]. A few problems, however, remained unsolved. Thus, following an earlier work [15], Kargiet [5] 
mainly focused attention on the ratio [S] /|E|y at both ends and concluded again that QSSA can be implemented in either case, 


including the intermediate region. Bajzeret [7], on the other hand, maintained that QSSA is valid, if at all, only for a short time- 
interval when [E], |S pis large and that the region of validity is considerably larger for large [S |, I|E|). The concentration profiles 
were calculated numerically and finding the steady state analytical expressions corresponding to the concentrations of substrate, 
enzyme- substrate complex, free enzyme and product. The purpose of this communication is to derive an approximate analytical 
expression for the concentrations in terms of rate constant using Homotopy perturbation method [16- 18]. 


2. MATHEMATICAL FORMULATION OF THE INITIAL VALUE PROBLEM 


The two-step Michalies - Menten model is symbolized by the following reaction scheme 


ky ky 
E+S a ES—> E+ P 
k 


af 
(1) 


Where [S](t), [E](t), LES](t) and P(t) denote, the free substrate, free enzyme, the enzyme - substrate inter-mediate complex and 
product respectively. The parameters k_;,k, and ky are positive constants for each reaction. 


By applying the law of mass action, which states that reaction rates are proportional to the concentrations of the reactants, the 
time evolution of scheme (1) can be determined from the solution of the system of following coupled nonlinear differential 
equations: [7] 


AS) <,_fes]-1 {es @ 
AE] (4. +e,)[es}-4,{6s] 3 
AES] _  fel[s]-(4 +4.) A 
alP] = ky [ES] 5) 


with the initial conditions are: 

At t=0: [S]=[S]p. [E]=[E]p. 
[Es]=0 and [P|=0 

Since the enzyme E is a catalyst, its total concentration must be a constant. This conservation law is readily obtained by adding Eqns. 

(3) and (4): 


(6) 
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ale] , ales] _ 


0 (7) 

dt dt 
The integration of this equation leads to: 

[E\(t)+LESI(t) =[Elo (8) 
Similarly from eqns. (2), (4) and (5) we get the following relations: 
als] , a[és]  alP]_, je 

dt dt dt 
The integration of this equation leads to: 
[S](t)+[P1(t)+ [ES 1(t) =.5 19 (10) 
With this system of ordinary’ differential equations reduce to only _ two, for|S and [ES], namely 
d\s 

als) = -k,((£o1-[Es)[s]-&[£s] (11) 

d\|ES 

ales) =k, ((Eo]—-[ES][S]-(k. +k. ES] (12) 


3. ANALYTICAL EXPRESSION OF CONCENTRATIONS USING HOMOTOPY PERTURBATION APPROACH 

The analytical solution of non-linear equation is of great importance due to its wide application in scientific research. The Homotopy 
perturbation method (ref. Appendix A) [16-18] is used to give the approximate analytical solution of non-linear reaction/diffusion 
Eqns. (11) and (12) are 


ky 
k, Eot 4 
k_y + ky se ' 
[s] (t) =Sp exp k2Ep enkEo! o(kath) , (13) 
k | Ey (ki +k) 


Bo (-hFo" ~e{ei+k)) 


k 


[ES |(t)= 
where 
k=k_j +k, —-kyEo. (14) 
We have derived Eqns. (8) and (10) using Eqns. (14) and(15) we get the approximate analytical solution respectively. 
ie ( -KEyt _ o-(kitk, Ir} 

k (15) 


[E]®)=Ep 


- —k, Egt+ 
k_jtk, re 


[P(+)=50 - So exp K2Ey (eh Bot lkath)t 
k | hE (k4+h) (16) 


[see (o£ talk tk,)t } 


k 


4. NUMERICAL SIMULATION 

The non-linear differential Eqns. (2-6) for the given initial conditions is being solved numerically. The function matt, in Matlab 
software is used to solve the ordinary differential equations given in Appendix B. The numerical solution is compared with analytical 
solution in Figs. 1-3. Satisfactory agreement is noted between the analytical and numerical results. 


5. DISCUSSION 

Eqns. (13)-(16) are the closed and simple analytical expression of concentrations of the free substrate, enzyme - substrate inter- 
mediate complex, free enzyme, and product respectively. Figure (1-3) describes the dimensionless unsteady-state concentrations 
versus time ft .The analytical expression of concentrations of the free enzyme[£](t), free substrate S(t), the enzyme - substrate 
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inter-mediate complex [ES|(t) and the product[P](t) have been plotted for various values of the parameters k_;,k,; and ky. From 


these figures it is inferred that the concentration of substrate [S](t) follows a first order exponential decay and it is always 
decreasing function from the initial value. The concentration of enzyme-substrate complex [ES(t) and product[P](r) initially 
increases and attains its steady state value at short intervals of time for all values of parameters. Free enzyme [E](t) increases slowly 


and reaches the steady state when time is very large. The time taken the steady value depends upon the values of parameters 
k_1,k, and ky . Our approximate analytical expressions of free substrate, enzyme substrate and free enzyme concentration are 


compared with simulation results in figs.1-3. A satisfactory agreement is noted. 


(a) 


[Slo = 1[E], =10 


[ES]o = [Php = 0 
ky kyky =0.1 


centrarion[S(t] [Ele [ESir] (Pir) 


Concentration|[S(r)] [E(e} [ES(e | [Pir)] 


Cone 
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(©) 


[s}, =U{E]. =10 
[Es], = [Ph =0 


zs, AE =10 


ConcenrratiorfS(r] [Er] [ES(r] [Pir] 
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Figure 1 Plot of concentrations [S()L. [E(¢)] [ES(t)] and [P(t)] versus time? for different 
values of the parameters k_},k and k3 . The curves are plotted using Eqns. (13 - 16), for some 
fixed values of the other parameters: 
) 0 and k_,,k>,k3=0.1 
(6) [S] =1, [El =10, [ES], =[P] =0 and k_,k..k3=1 
O and k_,,k,k;=10 


6. CONCLUSION 

The times dependent non-linear equations have been solving analytically. The obtained results have a good agreement with those 
obtained using numerical method. The primary result of this work is simple, straight forward and a new method of estimating the 
concentrations of substrate, product, enzyme- substrate complex and enzyme are derived for all possible values of parameters 
k_,, ky, and kg. The solution procedure can be easily extended to all kinds of system of coupled nonlinear equations with 


various complex boundary conditions in enzyme — substrate reaction diffusion processes. The analytical method is an extremely 
simple method and it is also promising method to solve other nonlinear equations. 
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Appendix A. The analytical solutions of non-linear Eqns. (11) and (12) 
In this appendix, we derive the solution of non-linear equation Eqns. (11) and (12) using new approach of Homotopy perturbation 
method. 


The Homotopy for the nonlinear equation (11) and (12) can be constructed as follows 


d\S 
a-p {I +ky(Eo 115 


(_- pes) —kLEolIS1+ (K+ tole + 
(A.1) dt 


i p{ +hy(Ey1-[Es)ISI-k es} =0 p{ eh ate _[es)isi—(k, +k] ~9 


(A.2) 


dt 


[Slo = L[E], =10 
[ESlo =[P], = 0 


Concentration{S(t ] (Ele [Es(r] [Pir | 


Time (r) 
(©) 
= 
B 
aa [Sh =1[E]. =10 
= [E5]} = [Ph =0 
a 
= ae Ey kp =10 
5 
s 
& 7 
o 
ae o1 0.2 03 04 0. 0.6 o7 08 os 1 
Time (r) 


Figure 2 Plot of concentrations [S(t)L, [E(¢)I LES (t)] and [P(t)] versus time? for different values of 
the parameters k_,,k2 and k3 .The curves are plotted using Eqns. (13 - 16), for some fixed values of 
the other parameters: 

0, [ES] =[P]p =0 and k_1,k2,k3=0.1 

(>) [S]o =1, [Elo =10, [ES]p =[P]p =0 and k_y,k2,k3=1 

(C)[S]p =1, [E]p =10, [ES]p =[P]p =0 and k_1,k2 kz =10 


Suppose the approximate solutions of Eqns. (A.1) and (A.2) have the form 
[s]= [slo + pls, + p’[s], +... (A3) [ES |= [es ]o + ples], 4 [es], +... (A4) 


Substituting Eqns. (A.3) and (A.4) into Eqn. (A.1) and Eqn. (A.2) respectively equate the terms with the identical powers of p, we 
obtain 


d 

pe: Ph ~ki [Eo ILs]o =9 a 
d 

ae dlesle —kLEg [slo + (ky + ko Jles]p =0 (A6) 
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The initial conditions are as follows: 


t=0: [s]y = [Sls and [es |p =0 (A.7) 
Solving the Eq. (A.5) using the boundary conditions Eq. (A.7), we get 
[s]p (4) = Spe" (A.8) 


Solving the Eq. (A.6), substitute Eqn. (A.8) in Eqn. (A.6) and using the boundary conditions Eq. (A.7), we get 


-kEyt — -(k.+k,)t 
eh atta ae) Ag 
where k=k_, +k, —k, Eq 


Substitute Eqns. (A.8) and (A.9) in Eqn. (11) we get 


[5] (.) =Sp exp k2Ey{ ek! olkitk,) 1) |(A-10) 
kiEy — (k. +ka) 
where k=k_, +k, —k, Eg 


Higher order iteration will be considered to improve the accuracy of the results. 


(b) 


= 
Ss 


© 
© 
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Figure 3 Plot of concentrations [S(t)], [E(¢)] L[ES(t)] and [P(t)] versus timef for different 
values of the parameters k_},k2 and k3 .The curves are plotted using Eqns. (13 - 16), for some 
fixed values of the other parameters: 

(a) [S]y =1, [E] =10, [ES] =[P], =0 and k_,k,k3=0.1 

(6) [S]p =1, [Ely =10, [ES], =[P]p =0 and k_,,k2.k3=1 

(c)[S]) =1, [E]) =10, [ES]) =[P]) =0 and k_,,k,k3=10 


www.discoveryjournals.org OPEN ACCESS 


disc@ very 


Paged 2 


ARTICLE 


Appendix B: Scilab/Matlab program for the numerical solution of nonlinear differential equation (5) 


function matt 
options= odeset('RelTol’,1e-6,'Stats','on’); 
%initial conditions 


x0 = [1;10;0;0]; 
tspan = [0,5]; 
tic 


[t,x] =ode45(@TestFunction,tspan,x0,options); 
toc 

figure 

holdon 

plot(t, x(:,1)) 

plot(t, x(:,2)) 

plot(t, x(:,3)) 

plot(t, x(:,4)) 

plot(t, x(:,5)) 
legend('x1','x2','x3','x4','x5') 
ylabel('x') 

xlabel('t') 

return 

function [dx_dt]= TestFunction(t,x) 
k1=1;k2=1;k3=1; 

dx_dt(1)= k3*x(3)-k1*x(1)*x(2); 
dx_dt(2) =(k3 +k2)*x(3)-k1*x(1)*x(2); 
dx_dt(3)=k1*x(1)*x(2)-(k3 +k2)*x(3); 
dx_dt(4)=k2*x(3); 


dx_dt = dx_dt’; 
return 
Appendix C 


Nomenclature 


Symbol Meaning Usual dimension 
[s] Concentration of free substrate moll 
[E| Concentration of free Enzyme Emolmin-! 
Concentration of free Enzyme — Substrate 
intermediate complex e4 
4 
[P| Concentration of Product mol/l @ 
[sh Initial concentration of free substrate moll 2 
[El Initial concentration of free Enzyme mol/l 
Initial concentration of free Enzyme — Substrate 
[ES] : . mol/l 
intermediate complex 
[Plo Initial concentration of Product mol/l 
k_4,ky,k2 Positive rate constant for each reaction mot! min! 
t Time factor min ae) 
+ 
85 
z 
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